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Abstract 

We calculate the new physics contributions to the neutral B'^ and B'^ meson mass splitting 
AMrf and AMg induced by the box diagrams involving the charged-Higgs bosons in the top 
quark two-Higgs doublet model (T2HDM). Using the precision data, we obtain the bounds on 
the parameter space of the T2HDM: (a) for fixed Mh = 400 GeV and 6 = [0°,60°], the upper 
bound on tan /? is tan (3 < 30 after the inclusion of major theoretical uncertainties; (b) for the 
case of tan /? < 20, a light charged Higgs boson with a mass around 300 GeV is allowed; and (c) 
the bounds on tan/3 and Mh are strongly correlated: a smaller (larger) tan/3 means a lighter 
(heavier) charged Higgs boson. 
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I. INTRODUCTION 



As a flavor changing neutral current process, B^—B^ mixing with q = d,s are generated 
at the loop-level and have been of fundamental importance in probing virtual effects from 
potential new physics beyond the standard model (SM). The B^ — B^ mixing is responsible 
for the small mass differences between the heavy and light mass eigenstates of neutral B 
mesons: 



AM, = M^o - M^o . 1 



The mass splitting AM^ has been measured with high precision 0, while the mea- 
surement of AMg is very difficult due to the rapid oscillation of Bg meson and has been 
reported by CDF and DO Collaboration 0, 0, 13 very recently. The world average for 
AM(i j3] and the first observation of AM^ from CDF 0| are the following 

AMd = 0.507 ± 0.005 ps-\ (2) 
AM, = 17.77 ± O.lO(stat) ±0.07(syst)ps-\ (3) 

which agree well with the standard model (SM) predictions or the results from global fit 
0. The perfect agreement between the SM prediction and the experimental measurements 
permit us to put strong constraints on the parameter spaces of various new physics models. 

In the SM, B^ — B^ mixing is dominated by the box diagrams with two internal t- 
quarks and W gauge bosons. In new physics models, the box diagrams with one or two 
W gauge bosons replaced by the new charged scalars or vector bosons and/or top quarks 
replaced by new fermions can also contribute to B^ — B^ mixing. Using the precision 
data, we study both B^ — B^ and B^ — 5° mixing in the top quark two-Higgs doublet 
model (T2IIDM) and try to find the constrains on the parameter space of this model. 

During the past years, B^ — B^ mixing has been studied extensively in the SM and 
various new physics models. The charged-Higgs boson contributions to B'^ — B^ mixing 
have been calculated at the leading order (LO) for a long time 0]. The next-to- leading 
order (NLO) quantum chromodynamics (QCD) correction to B^ — B^ mixing is firstly 
presented in Ref. jsl and the analytic formulae for the QCD renormalization group factors 
are given in Ref. [9|. In Ref. the authors studied the new physics effects in the 
conventional two-Higgs-doublet model (2HDM) of type I and II. In Refs. |il H [13 the 



authors have calculated the charged Higgs boson contributions to the mass splitting AM^^ 
and drawn the constraints on the parameters of the third type of 2HDM (the model III) 
at the LO or NLO level. Very recently, B^ — B^ mixing has been used to put constraints 
on various new physics models, for example, in Refs. jlj] after the release of the new date 
of AM,. In this paper, we will calculate the new physics contributions to AM, within 
the framework of T2HDM 3, 3 3 • comparing the theoretical predictions with the 



precision data, we draw the constraints on the free parameters of T2HDM. 

The organization of this paper is as follows. In sec. IHl we firstly give a brief review for 
the top quark two-Higgs doublet model, and then present the one-loop contributions to the 
mass splitting AM^ and AM, induced by the box diagrams involving the charged Higgs 
bosons in the T2HDM. Numerical results are presented in sec. IIIIl and the conclusions 
are included in the final section. 
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II. - B° MIXING IN THE T2HDM 

The new physics model considered here is the T2HDM proposed in Ref. 3] and studied 
for example in Refs. 3, 12, 3' which is also a special case of the 2HDM of type III (l9l |. 
This model is designed to accommodate the heaviness of the top quark by coupling it to 
a scalar doublet with large vacuum expectation value (VEV). All the other five quarks 
are coupled to another scalar doublet, whose VEV is much smaller. As a result, tan/5 is 
naturally large in this model. 

Let us now briefly recapitulate some important features of the T2HDM Consider 
the Yukawa Lagrangian of the form: 

Cy = -Ll<PiEIr - QL<PiFdR - qJ>^G1^^\r - qJ>2G1^^\r + H.c. (4) 

where Ql and are 3- vector of the left-handed quark and lepton doublets; (pi {i = 1, 2) 
are the two Higgs doublets with 0j = iT2(f)*; and E, F and G are the 3x3 matrices in the 
generation space and give masses respectively to the charged leptons, the down and up 
type quarks; l*-^-* = diag{l, 1, 0) and 1*^^^ = diag{0, 0, 1) are the two orthogonal projection 
operators onto the first two and the third families respectively. 

The Yukawa couplings involving the charged-Higgs bosons are of the form [il 

= -jM- — {-ULVMDdR[G+ - t&n(3H+] +URMuVdL[G+ - tan/?ff+] 

+URj:^VdL[t&nf3 + cot /3]H^ + h.c.}. (5) 

where G^ and denote the would-be Goldstone bosons and the physical charged Higgs 
bosons, respectively. Here Mjj and M^, are the diagonal up- and down-type mass matrices, 
V is the usual CKM matrix and S = MuUI^I^'^^Ur. is the unitary matrix which 
diagonalizes the right-handed up-type quarks and has the following form: 

1 




Ur = I sin0 COS0 I X I Vl-le^eP -ectr | • (6) 

e,ti Vl - \eM' 

where = rric/mt, and = |.^|e*'^ is a complex number of order unity. Inserting Eq. (jSj) 
into the definition of S yields 





m cejl^l' m,e,ievT^ra^ | • (7) 
m^eVl - le^eP mil - lectin 

Now we are ready to calculate the charged Higgs and would-be Goldstone bosons 
contributions to the mass splitting AMsgiq = d, s) in the T2HDM. 

The effective weak Hamiltonian for AB = 2 processes beyond the SM can be written 

as 



1=1 i=i 
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with 

Q2 = im) , Q3 = [fnbl] 

Qa = im) , Q5 = {qnbi) [q^hi) , (9) 

where q = s, d, corresponding to the operators of Bg and Bd system respectively, Pl,r = 
(1 =F 75)72, and a, [3 are color indexes. The tilde operators Qi {i = 1, 2, 3) correspond to 
the ones Qi {i = 1,2,3) with opposite chirality. 

In T2HDM, there are two CP-even scalars {H^, h^), one CP-odd scalar (A^), two 
charged Higgs bosons (i^^), and the Goldstone bosons {G^, G^). At one loop level only 
the charged scalars are relevant for the box diagrams contributing to the B^ — B^ mi xing 
amplitude. From Eq. (jSJ, we can rewrite the vertex couplings of = {H~^,G~^) [20| 
(where G^ is the would-be Goldstone boson) in a compact form 



Cnt = H^uaVai {af'PL + ai^'Pn) di + h.c. (10) 



where 



^AIl 



^^-1 tan/3, for 1=1 



^/2swMw I 1, for 1=2 




:ii) 



All _ e mj J tan/3, for 1=1 , , . 

V2swMw ■ I -1, for 1=2, 

with mA = {mu,mc,mt) and mi = {mdymsyrrib). The contributions of to the Wilson 
coefficients Ci of the relevant operators responsible for B^ — B^ mixing can be easily 
expressed in terms of the coefficients a;^^' and follows. 

The contributions to the Wilson Ci(/i) and C2(/i) induced by the box diagrams with 
one and one propagator can be written as ^: 

o2 



CiW = lyTL^ Tr-2 mAVfiA'a}^ ■ Doim^^m^,, M^, Mh+), 

A,A' ^ 

CM = ^ E E ^^^^^^^f^^'^ r^'^ • ^D,,{m% ml„ M^, <,), (13) 
where sw = sin6'vy {9w is the Weinberg angle), while the four-point integral functions Dq 



^ The contribution of to Ci{^) is already taken into account in the Inanii-Lim function So{xt) |2l|. 
Mass of the u quark are neglected. 
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TABLE I: The "magic numbers" appearing in the calculation of the Wilson coefficients in the 
process of 5° - mixmg. 





1 


2 3 


4 


5 




1 


2 


3 


4 


5 




0.286 


-0.692 0.787 


-1.143 0.143 








0.865 











(11) 

(22) 
(23) 
(32) 
(33) 
(44) 
(45) 
(54) 
(55) 


-0.017 














,(22) 





1.879 0.012 











-0.18 


-0.003 








,(23) 





-0.493 0.18 











-0.014 


0.008 








,(32) 





-0.044 0.035 











0.005 


-0.012 








^(33) 





0.011 0.54 











0.000 


0.028 








,(44) 








2.87 














-0.48 


0.005 


,(45) 








0.961 


-0.22 











-0.25 


-0.006 


,(54) 








0.09 














-0.013 


-0.016 


,(55) 








0.029 


0.863 











-0.007 


0.019 



and Dqq can be written as 

1 



L>o(a, 6, c, d) 



iir"^ (g^ — a)(g^ — — c)(g^ — d) 



a 1 r'^n ^ 1 

■ logi^J + 7 Tvn 77^ K log[ 



{b — a) (c — a) {d — a) d (a — h) ih — c){b — d) d 

log[^] , (14) 



Doo{a,b, c, d) 



(a — c) (6 — c) {d — c) d 
1 f d'^q (f 



4 j ivr^ (g^ — a)(g^ — — c)(g^ — d) 

a} a b^ b 

{b-a){c-a){d-a) ^ {a - b){b - c){b - d) 



c2 



log[3] • (15) 



(a — c){b — c) (d — c) d 

The contributions induced by the box diagrams with two Hj^ propagators can be 
written as ^: 

'J JAIl AJnJA'In A'Jl 
^AN - -^l^l^ ^2 «L «L «L 

A,A' l,n 

■Doo{m%ml„Ml^,Ml^y (16) 



^ In the sum over I and n in the expression for Ci(ii) the contribution of G^G^ is excluded since it has 
been taken into account in the function So{xt) (2l|. 
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n(,,\_ ^ ^AiVajVX'iVa'J JAIl^AJnJA'In^A'Jl 

'^iKN - ~2Z^Z^ iq:^ ^R ^R 

A,A' l,n 

■Doo (ml, ml, MI+,MI+^ , (17) 



CJlA - ^ W VXiVajVX'iVa'J JAIl AJn JA'In A'Jl 

A,A' l,n 



■Do (m%ml,,Ml+,Ml+^ 

(-^21/^; - ~2 2^ 2-^ iq:^ ^^a^u'^o-l «l % 

A,A' l,n 

■Do(m\,m\,,Ml^,Ml^y (19) 

2 

I \ "ST^ VXjVajVXi jVa' J iAIl AJn ]A' In A'Jl 

C^ifi) = -2^2^ m^Amu'^R a^' a[ a^' 

A,A' l,n 

■Do(m\,m\,,Ml^,Ml^y (20) 

2 

r(n\ - o ^AI^AjVX' I^A' J ]AIl AJn JA'In A'Jl 

'-si^j - ^2^2^ yq^^ ^ ^ ^ 

A,A' l,n 

■Doo (m\, m\,, M^+, M^+) . (21) 

These new physics contributions to the Wilson Coefficients are consistent with the ones as 
given in Refs. ji^l, the different sign and factor are due to the different definitions of the 
four-point integral functions. At onejoop level there are no contributions to the Wilson 
coefficients of the operators and Q^. 

Now, one needs to run the Wilson coefficients from the scale of new physics /i^ ~ My/ 
down to the low energy scale /if, ~ by using the QCDrenormalization group equations. 
For the evolution of these coefficients, we follow Ref. |2^, 



i=l s=l 



where we have set the new physics scale /ij = rrit and r] = as{fib)/as{mt). The magic 
numbers a,, bf'''^^ and c-^'*^ in Eq. (j22I) are listed in Table HI 

The off-diagonal element M12 in the 2x2 effective Hamiltonian causes the — 
mixing. The mass difference between the two mass eigenstates AM^ is described by 



(23) 



with 

M,? = (5j|K//(Ai? = 2)|i?°). (24) 
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In the SM, the mass sphtting AMg is calculated from the box diagrams of — 
mixing, dominated by t-quark exchange. At the NLO level, one fins that 23| 



(25) 



where Gi? is the Fermi constant, Mw the mass of the W boson, 5*0 (x^) the Inami-Lim 
function 2lj with Xt = m1/My^. The NLO short-distance QCD correction gives rjB = 
0.552, which is same for both B^ and B^ systems. The non-perturbative quantities Bb^ 
are the bag parameter and fs^ is the B^ meson decay constant ^24] . 

In terms of the bag-parameters, the matrix elements of the operators Qi and Qi are 
written as follows l2c 



\B,) 
\B,) 

\B,) 

\B,) 

\B,) 



lmBJlB['\^^l 



24 \mb{iJ,) + mq{jj,) 



24 \mb{n) + mg{fi) 



1 



mBJlBi''\^^) 



^BjlB^i\li) 



^Ba 



4 \mb{j^) + mq(/i) 
12 \mh{i2) + mg(/i) 



mBjlB^:\li)^ 
'niBjlA'\f^) 



(26) 



where Qi{fi) are the operators renormalised at the scale fi, Bi is the so-called bag factor. 
The matrix elements of Qis are the same as that of Q1-3. We use the same definition of 
B parameters B^''^ (/i) as in Ref. jl^ and find numerically that 



B['^\mi,) = 0.87(4) 
(m,) =0.82(3)(4), 
1.02(6)(9), 
1.16(3) 



h5 
-4) 



Bi'\m,) 
Bf\m,) 



B'^'imi,) = 1.91(4) 



+5 

-7' 
+22 
-7 1 



b[' 

-°2 

Bi' 
Bi' 



(m,) = 0.83(2)(4), 
(m,) = 1.03(4)(9), 
(m,) = 1.17(2)1^7, 
(mfe) = 1.94(3)lf . 



(27) 



III. NUMERICAL ANALYSIS 

In numerical calculations, we will use the following input parameters (all masses are 
in GeV) 

TUd = 5.4 X 10"^, TUs = 0.15, rrib = 4.6, 
rric = 1.4, mtijnt) = 165.9, m^^ = 5.279, m^^ = 5.367, 
v4 = 0.853, A = 0.225, p = 0.20 ± 0.09, r/ = 0.33 ± 0.05, (28) 

where A, A, p and f] are Wolfenstein parameters of the CKM mixing matrix. 
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TABLE II: The AMg{q = d, s) in the SM and the T2HDM for Mh = 300, tan /5 = 10, 30, 50 and 
5 = 0° (a), 30° (b) and 60° (c) within la range of the input hadronic parameters of JLQCD. 





SM 


T2HDM 


tan (3 = 10 


tan (3 = 30 


tan /? = 50 




t^in+o-oss 

U.01U_Q ]^33 


mBiitiiii 

(b)0.51llHi 


(a) 0.634lO:Ji5 

(b) 0.665l°:J2i 

(c) 0.569i°:}^^ 


(a) 1.630l0:l^ 

(b) 1.737t°i0? 

(c) 1.472l°:i0 


AMs 




(a) 17.15t|°9^ 

(b) 17.12l|06 

(c) 17.07i;°| 


(a) 20.09+|ff 

(b) 18.36l|26 

(c) 14.52l|63 


(a) 44.76l^:6i 

(b) 39.32l6j5 

(c) 26.06lt27 



For the hadronic parameters f^^Bsg, we use the values as given by the JLQCD col- 
laboration (26^ . 

fB,B'£UQCB = (0.2151°;°^^) GeV, 

/b.^'Ijlqcd = (0.245l°:°^J) GeV, (29) 

where the individual errors given in Ref. [i^ have been added in quadrature. 

In a previous paper fl^, we studied the new physics contributions to the B — > 
decay, and found strong constraints on the free parameters of the considered T2HDM: 

(i) A light charged Higgs boson with a mass less than 200 GeV is excluded. For 
fixed tan/3 = 30 and 6 = 0°, the lower limit on Mh is Mh > 300 GeV. 

(ii) The data of B Xg'y prefer a small angle 6: 6 < 44° for tan/3 = 30 and 
Mh = 400 GeV. 

Here we will consider these constraints in our choice for the free parameters of the T2HDM. 

The new physics (NP) contributions to B^ — B^ mixing are in general theoretically 
clean to interpret and have simple operate structure. To constrain deviations from the 
SM in these processes, we use the well measured physical observable AM^ as well as the 
first observation of AMg as given in Eqs. (|2l3|) and consider the effects of the hadronic 
uncertainty. The theoretical predictions about the mass difference AMg in the SM and 
T2HDM are listed in Table HIl for tan/3 = 10, 30, 50, M^^ = 300 GeV and 5 = 0°, 30° and 
60°. It is clear that the new physics contribution to B^ — B mixing is not sensitive to 
the parameter 6 when tan/3 < 20. 

In Fig. ^ we show the tan /3 dependence of AM^ and AMg in the SM and T2HDM with 

the central values of hadronic parameters. The band between two horizontal dot-dashed 

lines shows the measured values within 3a errors: 0.492 < AM^ < 0.522 (ps~^), and 

17.41 < AMg < 18.13 (ps~^). The solid horizontal lines show the central values of the 

SM predictions, which agree well with the data. The dominant theoretical error comes 

^ 1/2 

from the large uncertainty of hadronic parameter Jb^B^^ . The three short-dashed curves 
in Fig. ^represent the theoretical predictions in the T2HDM for Mh = 400GeV and for 
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FIG. 1: The plots of AMg vs tan (3 in the SM and T2HDM for Mh = 400 GeV, 5 = 0°, 30° and 
60°. 

5 = 0°, 30° and 60°, respectively. From this figure, the upper bound on tan/3 can be read 
off, 

tan/3< 25 (30) 

for Mh = 400 GeV and 5 = [0°,60°]. This bound is much stronger than the one ob- 
tained from the radiative decay B — >■ Xs'j After the inclusion of the effects of the 
uncertainties of the hadronic parameters, the upper bound on tan j3 will be changed into 

tan/5<30. (31) 

In Fig|2l we show the tan (3 dependence of AMg in the SM and T2HDM for 6 = 0° and 
for tan j3 = 10, 20, 30 and 40, respectively. Here the central values of the input parameters 
are used. The same as Fig. Q the band between two horizontal dot-dashed lines shows 
the measured values within 3a errors. The solid horizontal lines also show the central 
values of the SM predictions. The four short-dashed curves are the theoretical predictions 
in the T2HDM for fixed 6 = 0°, and for tan/5 = 10, 20, 30 and 40, respectively. It is easy 
to see that the data prefer a small tan f3, say tan (3 < 20, if one assumes the existence of 
a light charged Higgs boson with a mass around 300 GeV. The bounds on tan /? and Mh 
are indeed strongly correlated: a smaller (larger) tan (3 means a lighter (heavier) charged 
Higgs boson. 



IV. CONCLUSIONS 

In this paper, we have calculated the new physics contributions to the neutral B meson 
mass splitting AM^ and AMs induced by the one-loop box diagrams involving one or two 
charged Higgs boson propagators in the framework of T2HDM. 

By comparing the theoretical predictions with the precision data of — -B^(^-) mixing, 
strong constraints on the free parameters of T2HDM can be obtained. From the numerical 
results presented in last section, one can see that: 
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0.4 
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400 500 600 

MJGeV] 



400 500 600 

M„ [GeV] 



FIG. 2: The plots of AMg vs Mh in the SM and T2HDM, for 5 = 0° and for tan/3 = 10, 20, 30 
and 40, respectively. 

(i) For fixed Mh = 400 GeV and 5=[0°, 60°], the upper bound on tan /5 is 

tan/3< 30 (32) 
after the inclusion of major theoretical uncertainty. 

(ii) For a small tan /?, say tan (3 < 20, a light charged Higgs boson with a mass 
around 300 GeV is allowed. The bounds on tan/? and Mh are indeed strongly 
correlated: a smaller (larger) tan (3 means a lighter (heavier) charged Higgs boson. 
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